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Abstract. Long gamma-ray bursts are thought to be caused by a subset of exploding Wolf-Rayet stars. We argue that the 
circumstellar absorption lines in early supernova and in gamma-ray burst afterglow spectra may allow us to determine the main 
properties of the Wolf-Rayet star progenitors which can produce those two events. To demonstrate this, we first simulate the 
hydrodynamic evolution of the circumstellar medium around a 40 Mq star from the creation and evolution of a wind-blown, 
photo-ionized bubble around the star up to the time of the supernova explosion. Knowledge of density, temperature, and radial 
velocity of the circumstellar matter as function of space and time allows us to compute the column density in the line of sight to 
the centre of the nebula, as a function of radial velocity, angle, and time. While without radiative transfer modeling and without 
detailed knowledge of the spatial distribution of chemical elements we carmot produce spectra, our column density profiles 
indicate the possible number, strengths, widths, and velocities of absorption line components in supernova and gamma-ray 
burst afterglow spectra. Our example calculation shows four distinct line features during the Wolf-Rayet stage, at about 0, 50, 
150-700, and 2200 km s"', with only those of the lowest and highest velocity present at all times. The 150-700 km s"' feature 
decays rapidly as a function of time after the onset of the Wolf-Rayet stage. It consists of a variable number of components, and, 
especially in its evolved stage, depends strongly on the particular line of sight. A comparison with absorption lines detected in 
the afterglow of GRB 021004 suggests that the high velocity absorption component in GRB 021004 may be attributed to the 
free streaming Wolf-Rayet wind, which is consistent with the steep density drop indicated by the afterglow light curve. The 
presence of the intermediate velocity components implies that the duration of the Wolf-Rayet phase of the progenitor of GRB 
021004 was much smaller than the average Wolf-Rayet life time, which strongly constrains its progenitor evolution. 
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1. Introduction 

Planetary nebulae are the result of wind-wind interaction, with 
a fast wind emitted by a hot post-AGB star sweeping up a pre- 
viously emitted slow outflow (e.g., Villaver et al. I2002> . Stars 
above ^ 25...30Mo also undergo a red-blue evolution in the 
Hertzsprung-Russell diagram, moving from the red supergiant 
or Luminous Blue Variable (LBV) stage into the Wolf-Rayet 
phase (Meynet & Maeder l2000i) . As, in both cases, dense out- 
flows occur first at low and then at high speed, Wolf-Rayet stars 
are expected to produce circumstellar nebulae just like post- 
AGB stars (Garcia-Segura et al. '1996a' and 1996b'l, which are 
indeed observable in a large fraction of them (Miller & Chu 

Unlike AGB stars, massive stars end their evolution in a 
violent event triggered by the collapse of the massive iron core 
in the stellar interior While it cannot be excluded that some 
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of the Wolf-Rayet stars quietly collapse into a black hole, it is 
believed today that Type Ib/c supernovae and long gamma-ray 
bursts are both in fact exploding Wolf-Rayet stars. In Type Ib/c 
supernovae, photospheric hydrogen appears to be absent, while 
their location and their late time spectra are consistent with core 
collapse events (van Dvk ll992l Clocchiatti et al. 119961 Wang 
et al. ^OOSV Also hydrogen-rich Wolf-Rayet stars are potential 
supernova progenitors (Langer 1991 1. 

For gamma-ray bursts. Type Ic supernovae were found to 
dominate the late afterglow emission in two cases (Galama et 
al. ll998. Stanek et al. ,2003i . consistent with the predictions of 
the collapsar model for long gamma-ray bursts (Wooslev I 19931 
MacFadyen & Woosley'1999). Rigon et al. ( 2003) argue, based 
on possible GRB associations with the peculiar Type II super- 
novae SN 1997cy and SN 1999E, that some hydrogen-rich su- 
pernovae may also be related to gamma-ray bursts. 

In the supernova and the afterglow cases, the explosion cre- 
ates a powerful source of UV...IR photons, either in the super- 
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nova photosphere or in the jet interaction with its surroundings. 
While in gamma-ray bursts the photons produced by the inter- 
action of the explosively ejected material with the circumstellar 
matter — e.g., the gamma-ray burst afterglows — probe the in- 
nermost parts of the wind bubble (Chevalier et al. 20Q4), it is 
the absorption of photons in the line of sight to the central light 
source which traces the structure of the whole bubble and thus 
carries full information about the evolution of the progenitor 
star 

With the optical and near-infrared spectroscopy of the af- 
terglow of GRB 021004 starting only 0.6 days after the burst, 
it has become possible for the first time to use absorption line 
systems to analyze the entire circumstellar bubble surround- 
ing a gamma-ray burst, as proposed by Schaefer et al. ( I2003> . 
Mirabal et al. (2003 1, Fiore et al. (2004 1 and Starting et al. 
( 12005 > . With continuously refined GRB and afterglow observ- 
ing techniques, more and better afterglow spectra might pro- 
vide more such cases soon. 

Since Type II and Type Ib/c supernovae also have mas- 
sive star progenitors, one may also expect narrow circumstellar 
lines in their spectra. Indeed, Dopita et al. (1984) attributed 
narrow P Cygni profiles seen in Ha and H<5 in the spectrum of 
the Type IIL SN 1984E to a Wolf-Rayet wind with a velocity of 
about 3000 km s^' . Bowen at al. (^2000") and Fassia et al. j2UoTl l 
note the presence of a broad absorption feature at ~ 350 km s"' 
in the spectrum of the Type IIL SN 1998S, which, as they ar- 
gue, may be caused by a moving shell. Dedicated searches for 
circumstellar absorption lines in very early supernova spectra 
might again provide more and better cases. 

It is the main aim of this paper to demonstrate that circum- 
stellar lines in supernovae and GRB afterglows can be used 
not only to identify the Wolf-Rayet nature of their progenitor 
stars, but also to strongly constrain some of its main properties. 
Here, we present a 2D-hydrodynamic model of the evolution of 
the circumstellar medium surrounding a massive star, improv- 
ing upon the models produced by Garcia-Segura et al. (I I996al 
[l996b). Using this model, we then calculate the column density 
of the circumstellar medium as a function of the radial veloc- 
ity throughout the evolution, which allows a comparison with 
observed circumstellar spectral lines. 

This article is structured as follows: In Sect. |2 we give an 
overview of the method we use for our numerical simulation 
of the circumstellar medium. In Sects. |3and|3 we present the 
result of these simulations. In Sects. |5l and |6l we explain the 
method to compute the column density as function of radial 
velocity. Finally, in Sect.Qwe discuss the result of our simula- 
tions and in Sect.|S]we summarize our main conclusions. 

2. Numerical technique and assumptions 

At solar metallicity, stars in the mass range of about 25 Mq to 
40 Mq are thought to develop into red supergiants and there- 
after to become Wolf-Rayet stars. This means that three inter- 
actions take place in the circumstellar medium (CSM): 

1 . An interaction between the fast, low density main-sequence 
wind and the interstellar medium. 



2. The slow, high density Red Super Giant wind then hits the 
hot bubble created by the main-sequence wind. 

3. Finally, the massive, high velocity Wolf-Rayet wind 
sweeps up the previously created structures. 

For the hydrodynamical simulations, we used the ZEUS 3D 
code by Stone and Norman (1992 1. We simulated the main 
sequence and early red supergiant (RSG) phase only in ID. 
The Wolf-Rayet wind interaction was computed in 2D, so we 
took the end result of the ID simulation and map this onto a 
2D grid, as described in Garcia-Segura, Langer & Mac-Low 
(jT996Fl. However, since we followed the evolution of the cir- 
cumstellar medium up to the supernova stage, we calculated 
the whole circumstellar bubble in 2D, instead of cutting out the 
inner part. For simplicity, we divided the stellar evolution into 
three stages: main sequence, red supergiant and Wolf-Rayet. 
As input model we used the 40 M© star model with metallic- 
ity 0.02 as calculated by Schaller et al. d 19921 1. The resulting 
parameters for the stellar wind and photon count are given in 
Tabled 

Mass-loss rate, wind velocity and the number of ionizing 
photons were taken as average over these three stages. The 
mass-loss rate follows directly from the evolutionary model, 
as does the photon count, for which we used a black body 
approximation. The wind velocity was chosen so that the ki- 
netic energy output of the star during each stage is reproduced. 
We chose to use the kinetic energy rather than the momentum, 
since wind driven shells are energy driven for the largest part 
of their existence (Weaver et al. ll977l l. The density of the inter- 
stellar medium was set at 20"^^'^ g cm"^. The effect of photo- 
ionization was included in the simulation by calculating the 
Stromgren radius along each radial gridline and correcting the 
temperature and mean particle weight within this radius as de- 
scribed by Garcia-Segura & Franco (1996 1 and Garcia-Segura 
et al. ( I1999I I. A similar calculation for a 60 Mq star was re- 
cently performed completely in 2D by Freyer et al. ( <2003> . 

3. The CSIVI during main sequence and RSG phase 

From the beginning of the main sequence, the H II region cre- 
ated by radiation from a massive star pushes a shell into the 
interstellar gas. This shell consists of neutral hydrogen and is 
produced by a D-type front; i.e., the gas in the shell enters the 
ionization front subsonically. At the same time the kinetic en- 
ergy of the wind is converted into thermal energy through the 
collision with the surrounding matter. This creates a hot bub- 
ble which pushes a shell into the H II region. Eventually, the 
H II region and the wind-blown bubble reach pressure equilib- 
rium. Together they drive a shell into the interstellar medium. 
The shell between the wind blown bubble and the H II region 
disappears, though the density discontinuity remains. 

This sequence of events depends on the pressure balance 
between wind driven bubble and HII region. If the pressure 
in the wind-driven bubble becomes so high that it drives the 
shell into the H II region with supersonic speed, then the two 
regions will not find a pressure balance. In that case the wind- 
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Phase 


End of phase [yr] 


duration [10^ yr] 


m [ Mo yr '] 


Vwind [kms '] 


'^photon [s ] 


Main Sequence 


4.309x10* 


43.1 


9.1x10-^ 


890 


4.62x10-*' 


Red supergiant 


4.508x10* 


1.99 


8.3x10-^ 


15 


3.00x10*' 


Wolf-Rayet 


4.786x10* 


2.78 


4.1x10^^ 


2160 


3.86x10^*' 




R[pc] 



Fig. 1. Structure of the circumstellar bubble around our 40 
star at the end of the red supergiant phase. This figure gives 
density (continuous line) and thermal pressure (dotted line) as 
a function of the radius. From left to right we have: the freely 
expanding RSG wind, the wind termination shock with the thin 
RSG shell, the hot bubble, the energy driven shell, the former 
H II region, the 'old' shell (which was driven into the interstel- 
lar medium during the main sequence phase and is now dissi- 
pating, since the photo-ionization no longer extends beyond the 
wind bubble), and the interstellar medium. 



driven shell will sweep up the H II region. The end result will 
be a single bubble without density discontinuity, driving a shell 
into the interstellar medium. The density profile will resemble 
that shown in Garcia-Segura et al. ( i996b ), although the photo- 
ionization will cause instabilities in the shells (Garcia-Segura 
& Franco [19961 and Garcia-Segura et al. 1999 1. The deciding 
factors here are the ram pressure of the stellar wind, the number 
of ionizing photons, and the density of the interstellar medium, 
which together determine the thermal pressure in the bubble 
and the radius of the suiTounding HII region. Since the ram 
pressure of the wind depends on both wind velocity and mass 
loss rate, which in turn depend on such factors as stellar mass 
and metallicity, the ram pressure can vary considerably. The 
same is true for the number of high energy photons and the 
density of the interstellar medium. These two factors determine 
the radius of the H II region. If this radius is large enough, the 
wind-driven shell will have slowed down to subsonic speeds 
before it reaches the photo-ionization driven shell. In that case 
pressure balance wiU be reached after all, and a bubble with 



internal density discontinuity will be created. In this paper we 
shall continue to use the model with a two-part bubble, to show 
the widest range of possible phenomena in the circumstellar 
medium. 

During the red supergiant (RSG) phase, the number of ion- 
izing photons is low, so the HII region disappears. The wind- 
blown bubble keeps its high pressure and pushes a new shell 
into the suiTounding medium. A third shell forms at the loca- 
tion of the termination shock due to the lower ram pressure 
and much higher density of the stellar wind. At the end of the 
red supergiant phase, the circumstellar medium is built up as 
follows (see also Fig: 0: closest to the star is the freely ex- 
panding RSG wind, which ends at about 3 pc in the wind ter- 
mination shock. The shock itself is marked by a thin shell of 
shocked RSG wind material. Next comes the hot isotropic bub- 
ble, which in turn pushes a shell (r ~ 21 pc. in Fig.^ into the 
former H II region. Finally, at about 37 pc, there is the old shell, 
which was pushed into the interstellar medium during the main 
sequence phase, but is now slowly dissipating, since it is no 
longer supported by thermal pressure. The velocities of these 
three shells is low. In fact, the old shell at 37 pc is standing 
still, since there is no longer any force to drive it outward. The 
shell at 21 pc moves outward slowly (v ^ 10. .20 kms"'), driven 
by the thermal pressure in the hot bubble. The movement of 
the shell at the wind termination shock varies through the RSG 
phase. During the early part of this phase, the wind termination 
shock moves inward, due to the decreased ram-pressure of the 
stellar wind. The shell is formed during this phase. The ther- 
mal pressure in the hot bubble decreases over time, causing the 
wind termination shock (and the shell at that location) to move 
outward with a velocity of about 10 kms"'. For stars with a 
mass of 25 Mq or less, this is the end of the evolution, as was 
shown in van Marie, et al. ( I2004> . 

4. The CSM during the Wolf-Rayet phase 

The period that is of interest to us is the Wolf-Rayet phase, 
which can be subdivided roughly into three parts, based on the 
evolution of the circumstellar medium rather than the evolution 
of the star itself. During the first part (WRl), the ram pressure 
of the wind rises dramatically as the wind velocity increases 
(see Table [TJ. This creates a new shell, which moves rapidly 
through the free-streaming RSG wind and collides with the 
RSG shell at the wind termination shock (see Sect. O. Both 
shells are broken up during the collision, and their shattered 
remains coast into the surrounding bubble. This process was 
described in great detail by Garcia-Segura et al. ( 1996bb. The 
collision can be seen in Fig.|3 The sudden increase in the num- 
ber of ionizing photons causes most of the former H II region 
(see Fig.Q to re-ionize. 
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Fig. 2. The logarithm of the density [g cm '] of the circumstel- 
lar medium close to the star, just before (upper part, t - 4.514 
Myr) and just after (lower part, t = 4.522 Myr) the collision be- 
tween the RSG shell and the Wolf-Rayet shell. As can be seen, 
the Wolf-Rayet shell is already fragmented due to Rayleigh- 
Taylor instabilities before the collision. After the collision, the 
remnants of the two shells move outward together. 




20 30 

Gquatorial plana [pc] 



Fig. 3. Circumstellar bubble around a 40 Mq star during the 
second part of the Wolf-Rayet stage (WR2). This figure shows 
the density of the circumstellar medium after 4.566 Myr. 
Moving outward from the star, we have the freely expanding 
wind, the wind termination shock, the hot bubble, which drives 
a shell into the surrounding H II region, and finally the old shell. 
The white line shows the Stromgren radius. The blobs of mat- 
ter inside the hot bubble are the remnants of the RSG shell and 
Wolf-Rayet shell after collision (cf. Fig.O. They are dissipat- 
ing into the surrounding medium. 



Fig. 4. Similar to Fig.|3] but at at 4.598 Myr. The Stromgren 
radius is receding as the Wolf-Rayet wind adds more matter to 
the cuxumstellar medium. 
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Fig. 5. Similar to Figs. |3] and |4] but during the third part of 
the Wolf-Rayet stage (WR3; t = 4.630 Myr). The Stromgren 
radius is now at the same distance as the wind driven shell. The 
remnants of the RSG shell and Wolf-Rayet shell have almost 
completely disappeared. 



The second part of the Wolf-Rayet phase (WR2) is a pe- 
riod of stabilization. The isobarity of the wind-blown bubble 
was severely disrupted during WRl, and it takes some time to 
re-establish a constant pressure, since the remnants of the col- 
lided shells move into the bubble at velocities close to the speed 
of sound. During WR2 these remnants travel outward into the 
bubble, where they eventually dissipate. The high density of 
the Wolf-Rayet wind makes it difficult for the ionizing photons 
to penetrate deeply into the circumstellar bubble, causing the 
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Stromgren radius to move inward. The Stromgren radius will 
eventually end up at the same radius as the shell that is still be- 
ing pushed outward by the wind-blown bubble (see Fig.l^and 
13 . As the kinetic energy of the Wolf-Rayet wind is very high, 
the thermal energy of the wind-blown bubble will increase, 
causing the shell that it pushes to accelerate to a velocity of 
about 50 kms"'. 

During the third part of the Wolf-Rayet phase (WR3), there 
is only one hot bubble maintained by both the stellar wind 
and the photo-ionization. This bubble pushes a shell into the 
surrounding medium (former HII region). Photo-ionization 
causes the shells to fragment (Garcia-Segura et al. 1999 1, i.e. 
the radiation will penetrate through the moving shell at some 
location and start ionizing the material beyond. The result is an 
extremely fractured Stromgren radius (see Fig.|5|i. As the pres- 
sure in the wind-blown bubble decreases and the moving shell 
sweeps up more material, it will decelerate. It reaches the orig- 
inal main sequence shell (see Fig.^ almost at the same time 
as the central star explodes as a supernova. Note that not every 
Wolf-Rayet star necessarily goes through all three phases. It 
depends on the total time the star spends as a Wolf-Rayet star, 
which varies greatly depending on the metallicity and initial 
mass of the star. 

5. CSM column densities 

A gamma-ray burst afterglow can be described as a point 
source in the center of our grid. In order to compare the column 
densities, we take a single radial line in this grid. Moving out- 
ward from the center and at each grid point, we take the local 
radial velocity and the column density (density x radial length 
of the grid cell). This gives us the column density as a func- 
tion of radial velocity, which is equivalent to the absorption as 
a function of the blue-shift. Typically, the column density as a 
function of the velocity will look like Fig.|6l Winds give nar- 
row lines, while shells produce a much broader profile. Since 
the circumstellar medium shows considerable variations, we do 
this for all radial lines in the grid and take the average, in order 
to make sure that all observable features are accounted for. 

Temperatures in the circumstellar medium can rise as high 
as r ^ 10^ K. Therefore, it is necessary to spread the column 
density given in each grid point over a velocity range, dictated 
by the Maxwell-Boltzmann distribution for the thermal veloc- 
ity in a single direction: 

— 

\27TkT I 



e 



-mv,^l2kT 



(1) 



with m the particle mass, k the Boltzmann constant, and T the 
local temperature. The probability that a particle has a velocity 
Dr is given by P(!;,)d!;r- In order to calculate the column den- 
sity as a function of velocity, we first calculate the Maxwell- 
Boltzmann distribution for the local temperature. Multiplying 
the distribution function by the mass density gives us the mass 
density of particles moving within a certain velocity interval. 
We use a velocity interval Av^ of 1 kms"', so the mass of gas 
particles with thermal velocities in the interval Ai>r is: 

pP(t;,-)A.,=pl^J e-'™''/2*^A.„ (2) 




1 



Velocity 



Fig. 6. Schematically, the column density profile produced by 
a typical wind- wind interaction. Here, a slow wind (1) is over- 
taken by a much faster wind (2). If the collision is supersonic, 
a shell (S) is created. Winds give narrow absorption lines, since 
all particles have about the same velocity; only thermal broad- 
ening may somewhat increase the width of the line. Shells, on 
the other hand, contain material that is slowing down, so the 
material will be spread out over a larger velocity range, creat- 
ing a much broader line. Wind-wind interactions tend to last 
only a short period of time, but a similar situation occurs when 
the material in region 1 is stationary. 



with Ay, equal to 1 kms~'. In the co-moving frame, the result 
is a Gaussian curve with mean velocity zero; i.e., the amount of 
gas moving away from us equals the amount of gas moving to- 
ward us. We then shift this entire curve until the mean velocity 
equals the radial velocity of the gas, as found in the hydrody- 
namical simulation. Multiplying this by the length of a grid cell 
(dr) gives us the column density of mass in a velocity interval 
as a function of the velocity for one grid cell. We repeat this 
process for all the grid cells along a radial line and integrate 
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them for each velocity interval separately. The quantity we end 
up with is dc(u, Av): 

d,(v,Au)^ I p(r:)P(ur)AurdT, (3) 

Ji=0 

with R the total radial size of the grid. This is the column den- 
sity resulting from all particles along a radial line moving with 
velocities within 0.5 kms ' of u^. 

5.1. Angle-averaged column densities 

Figure shows the angle-averaged CSM column density as 
function of the radial velocity over the entire Wolf-Rayet pe- 
riod; i.e., in order to ensure that all features are contributing, we 
repeat the process described above for all 200 radial grid lines 
and take the average. To obtain the total column density within 
a given velocity interval, we need to specify the minimum and 
maximum velocity and then add up the column densities of the 
individual 1 kms"' intervals between those two. 

The chosen time interval in Fig. starts at the end of the 
red supergiant phase and ends at the pre-supernova stage of 
the underlying stellar model. While the central star is a red su- 
pergiant, there is only one clear line feature, at a radial veloc- 
ity of nearly zero. There are no high speed features, since the 
RSG wind, the shells, and the gas in the hot bubble, all move at 
low velocities. As the Wolf-Rayet wind starts, two new features 
appear: 1) the Wolf-Rayet wind itself, traveling at about 2200 
kms"\ and 2) the region where the Wolf-Rayet wind collides 
with the much slower RSG wind (i.e. the Wolf-Rayet shell). 

During the initial Wolf-Rayet phase (WRl in Sect.0} three 
absorption features are produced: the zero velocity feature, 
which was already visible during the RSG phase, a narrow 
strong line at high velocity (2200 km s ') caused by the Wolf- 
Rayet wind, and a feature at intermediate velocity, which is 
first produced by the swept-up red supergiant wind shell and 
later by its collision with the RSG shell (cf. Fig. 2). Initially, 
the intermediate velocity feature extends over a large velocity 
range (100. ..600 kms '). This is due to both the high temper- 
ature (which causes the particle velocities to be spread over a 
large velocity range) and the fact that gas moving at a veloc- 
ity of 2200 km s ' is slowed down abruptly. However, this fea- 
ture becomes narrower over time. The bulk material of the shell 
moves outward at a velocity of about 200 kms"' and eventu- 
ally hits the RSG shell (see Fig. After the shell collision, 
the largest fragments maintain their original speed, which is 
between 150 and 300 kms"'. In some places shocked Wolf- 
Rayet wind material breaks through the fragmented shells to 
form smaller clumps, which can have velocities of up to ca. 
700 kms"'. This marks the beginning of phase WR2. 

During this phase the number of visible features is variable. 
There is still the intermediate velocity feature caused by the 
remnants of the RSG and Wolf-Rayet shell. However, since this 
feature is then caused by shell fragments which cover only a 
fraction of the possible viewing angles (cf. Figs. 0] and |5}, it 
may or may not show up in an observation (see below). As the 
shell fragments start to dissipate into the surrounding medium, 
the intermediate velocity feature disappears. 



At about the same time, a new feature appears close to the 
zero velocity line. This feature, which shows a velocity of about 
50 kms"', is the result of the wind driven shell which moves 
outward as the pressure in the wind blown bubble increases (see 
Sect. 1^. Since this shell is not fractured, the line will always 
be present, although it may be difficult to identify due to its 
proximity to the zero velocity feature. 

The period during which the intermediate velocity feature 
exists and the 50 kms"' feature is formed is shown in Fig. 
m which is a magnified part of Fig.0 It shows that the inter- 
mediate velocity feature is concentrated between 150 and 300 
kms"', with subcomponents of up to ca. 700 kms"' present 
during the first 40 000 yr of the Wolf-Rayet stage. Figure |9] 
shows vertical cuts of Fig. 0for the times which correspond 
to the CSM density plots in Figs. |3] to |5] During the third 
WR phase, WR3, the 50 kms"' feature disappears as the 
shell moves further outward, sweeping up material and slow- 
ing down. In the end only two lines are visible, one at zero 
velocity and one at the velocity of the Wolf-Rayet wind. 

Any ion which is present throughout the circumstellar 
medium should show up in a set of absorption features as de- 
scribed here. The relative height of the lines, or depth of the 
absorption, is more difficult to estimate as it depends on the 
percentage of the particle density of the gas that the ion takes 
up, which depends both on the chemical composition of the gas 
and on the degree of ionization (cf. Sect. 17. 3> . 

5.2. Line of sight dependence 

The inclination of the faster parts in the intermediate velocity 
feature in Figs. 0and|8] resp. — i.e. the fact that these fea- 
tures do not appear as vertical structures, but are somewhat 
inclined — may lead us to expect multiple absorption com- 
ponents at certain times. However, since the displayed column 
densities are averaged over all angles, the presence of multiple 
components of the intermediate velocity feature in the angle- 
averaged representation does not warrant visibility along a sin- 
gle line of sight (see below). On the other hand, even if the 
angle-averaged representation does not indicate multiple com- 
ponents for a given time, they might well exist at this time for 
specific lines of sight. We therefore provide a discussion here 
of the line of sight dependence of the various features. The zero 
velocity line, the Wolf-Rayet wind line, and the line caused by 
the intact shell (at 50 kms"') can be observed independent of 
the angle. 

The intermediate velocity feature caused by the fragmented 
Wolf-Rayet and RSG shells is only visible if a major fragment 
is in the line of sight. Its probability is strongly dependent on 
time, as can be seen in Fig. ^| In this figure we have taken 
the average column density of the gas moving at velocities be- 
tween 150 and 300 kms"', between 400 and 450 kms"', and 
between 450 and 600 kms"'. The first interval contains most 
of the shell fragments; the last contains the clumps with mod- 
erately high velocity (cf. Fig.|8}. If both these intervals contain 
a high column density while the intermediate interval does not, 
they may show up as separate absorption components in an ob- 
served spectrum. 
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As Fig. ^1 shows, the intermediate velocity feature is vis- 
ible at nearly any angle throughout the entire velocity interval 
before the Wolf-Rayet shell collides with the RSG shell (t ^ 
4.525 X 10^ years). After the collision, the shell remnants are 
broken up into independent, confined fragments. During this 
period the feature can only be observed if such a fragment hap- 
pens to be in the line of sight. 

The probability of such an occurrence is estimated in Fig. 
rm This figure shows the percentage of radial gridlines, along 
which the column density at 200 km s"' is more than five times 
as high as the column density at 100 kms"'. As can be seen, 
the chance of observing the line is highest before the shell col- 
lision; afterwards, it drops rapidly. The increased probability of 
observing the intermediate velocity feature at about 4.57 mil- 
lion years (Fig. II 1> is due to dissipation of the shell fragments. 
These fragments reach the edge of the hot bubble, at which time 
they are unable to continue their radial movement and start to 
spread out in an angular direction. However, by then the col- 
umn density has dropped appreciably, and the velocity tends to 
zero. 

Figure^Jclearly shows the complexity of the intermediate 
velocity feature in the angle and velocity space. While based 
on a 2D simulation with limited resolution, quantitative state- 
ments are difficult, it appears clearly from this figure that may 
expect to observe multiple velocity components in the interval 
of 100 to 700 km s"' for certain times and viewing angles. It is 
obvious that the clumps have considerable angular velocities, 
which means that they may drift in front of each other during 
their movements in the hot bubble. 

Figure ^] shows the column density as a function of the 
velocity along a single line of sight, at a time shortly after the 
collision of the WR shell and the RSG shell. Time and angle 
were chosen from inspection of Fig. [TOI such that the lower 
(150...300kms-^) and higher (450.. .600 km s"') velocity in- 
tervals of the intermediate velocity feature are populated, but 
the middle one (400.. .450 kms"') is not (pink color in Fig.llO>. 
And indeed. Fig. suggests that an observer looking at the 
specified time and angle might see two separate components 
of the intermediate velocity feature. However, it is also clear 
from Fig.Elthat the number of observed separate components 
may depend sensitively on the spectral resolution of the obser- 
vation. It should be emphasized that compared to Figs. ID 
and|9] which show averages over all angles, the structure of 
the intermediate velocity feature is much more complex when 
viewed along a single line of sight. 

6. Chemical composition of the CSIVI 

The column densities calculated in Sect. |5l are based on the 
assumption that a given ion is present throughout the circum- 
stellar medium. This, of course, depends on both the chemical 
composition and the temperature of the material. The composi- 
tion of the freely expanding Wolf-Rayet wind, which produces 
the 2200 kms"' feature, changes over time as the star evolves 
from the WN into the WC phase (Maeder & Meynet 1994). The 
composition of the wind will depend on the initial metallicity 
of the Wolf-Rayet star. During the WN phase, the wind will 
consist mostly of helium with some heavier ions. As the star 



becomes a WC star, carbon (produced by He burning) appears 
at the surface, becoming the dominant component of the wind. 

The wind blown bubble has a mixed composition, since it 
contains all the material that the star has ejected during its evo- 
lution. The remnants of the RSG and Wolf-Rayet shell, which 
cause the intermediate velocity line, have the composition of 
the RSG and early Wolf-Rayet phase, plus some material of 
the hot bubble that they have swept up as they travel outward. 
The main source of heavy elements here is the red supergiant 
wind. The energy driven shell on the edge of the bubble, the 
H II region, and the old main sequence shell, all have the com- 
position of the interstellar medium. 

The degree of ionization of the circumstellar gas depends 
on both temperature and photo-ionization. Ignoring photo- 
ionization, the temperature of the Wolf-Rayet wind would be 
comparatively low. The wind leaves the star with a tempera- 
ture of about 10^ K but would cool down rapidly. Ionization 
in this area comes from the radiation of the star. The temper- 
ature of wind-blown bubbles depends on the kinetic energy of 
the wind and density of the medium into which the bubble ex- 
pands. Before the Wolf-Rayet shell collides with the RSG shell, 
temperatures in the bubble that drives the Wolf-Rayet shell can 
rise as high as 10^ K. After the collision, the remnants move 
into the bubble created by the main sequence wind, where tem- 
peratures are lower. As the bubble becomes isotropic, the tem- 
peratures in the high density areas will be lower, but the temper- 
ature of the bubble as a whole increases due to the high kinetic 
energy of the Wolf-Rayet wind. Throughout the bubble, tem- 
peratures between 10^ K and 10^ K can be found. The energy 
driven shell and the area beyond form an H II region, which has 
a temperature of about 10"* K. The old shell and the interstellar 
medium are cold. This is the situation before the central star 
collapses. However, we should keep in mind that the emission 
from the gamma-ray burst afterglow passes through a medium 
through which the high energy photons of the gamma-ray burst 
itself have already passed, which is not taken into account in 
our model. 

7. Implications for GRB afterglow spectra 

Figures to ^2 show that circumstellar absorption features, to 
be expected in a gamma-ray burst afterglow spectrum, strongly 
depend on the duration of the Wolf-Rayet phase of the progeni- 
tor star. This may enable us to draw several conclusions regard- 
ing a gamma-ray burst progenitor and its surrounding medium. 

7. 1 . The progenitor star 

Absorption lines at zero velocity and the Wolf-Rayet wind 
velocity might be observed throughout the Wolf-Rayet phase. 
However, the absorption line at intermediate velocity is only 
visible for a short period of time (approx. 80 000 yr). The exact 
speed at which this line appears depends on the relative mass- 
loss rates and wind velocities during the Wolf-Rayet phase 
and the preceding stage. This line, caused by the Wolf-Rayet 
wind-driven shell, reaches its final velocity while the shell is 
expanding into the freely expanding RSG wind. After it hits 
the RSG shell at the wind termination shock, the high thermal 
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Fig. 7. Column density dc(ii, Av) with Au = I km s"' of the circumstellar medium around a Wolf-Rayet star as a function of radial 
velocity and time, averaged over 200 radial grid lines. On the horizontal axis: the time since the birth of the star in years. On the 
vertical axis: the radial velocity in km s ' . The plot starts at the end of the red supergiant phase, when there is only matter moving 
at low velocity. As the Wolf-Rayet phase sets in a narrow feature at the Wolf-Rayet, wind velocity appears as does a much broader 
feature at intermediate velocity, caused by the collision between Wolf-Rayet and RSG wind. This feature disappears after about 
80000 years, when the remnant of the Wolf-Rayet shell dissipates into the circumstellar bubble. The low velocity line splits at 
about the same time into two separate lines close together, as the increased thermal pressure of the wind bubble drives a shell 
into the surrounding medium. 



pressure in the hot bubble may slow down the shell fragments, 
although this does not happen with this particular model. 

It is possible to use an analytical estimate of the shell ve- 
locity during its earlier phase, in order to find the velocity 
at which the absorption line will appear. The velocity of an 
energy-driven shell, expanding into a previously emitted stellar 
wind, is given by Kwok t,2000j as: 



I 3mi 



(4) 



with Vs the velocity of the shell, Vi and ihi the velocity and 
mass-loss rate of the old wind, and V2 and ma the velocity and 



mass-loss rate of the new wind. With typical velocities and 
mass-loss rates for RSG and Wolf-Rayet winds, this yields a 
shell velocity on the order of 200 km s \ which is indeed what 
we found. However, both wind velocity and mass-loss rate de- 
pend on the metallicity and the initial mass of the star, and they 
can vary even for one star, if the stellar wind is not spherically 
symmetric. Since rapidly rotating stars have aspherical wind 
distributions and gamma-ray burst progenitors are thought to 
be rapidly rotating stars, this is a rational scenario. For heav- 
ier stars, which pass through a Luminous Blue Variable stage 
rather than a red supergiant stage, the values would be different 
as well (Garcia-Segura et al. ll996al . Whether the line will be 
visible at all depends ,of course, on the exact moment at which 
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Fig. 8. Blow-up of the low velocity part of Fig.0for the period during which the intermediate velocity feature is visible. This 
feature starts at about 4.508 Myr and is the result of the shell swept up by the Wolf-Rayet wind. Initially, the column density 
between the zero velocity feature and the intermediate velocity feature is quite high, as slow-moving RSG wind material is 
swept up by the shell and accelerated. After the collision between the Wolf-Rayet shell and the RSG shell (t - 4.52 Myr), this 
situation no longer applies. The fragments of the shells no longer sweep up a significant amount of mass. Therefore, no matter 
is accelerated and the column density between the zero velocity feature and the intermediate velocity feature decreases. The 
intermediate velocity feature is strongest between 150 and 300 kms"', but fainter components occur at velocities up to 700 
kms At the end of this period, the zero velocity feature splits up into two separate parts. The new feature, showing material 
moving at about 50 km s is created by the wind-bubble-driven shell, which moves out into the H II region. 



the star collapses. If the line is visible, it would indicate a short 
Wolf-Rayet phase, i.e. no more than about 50000 to 80000 
years. 

It is also possible that the freely expanding wind area 
around the star is much larger. Possible explanations for this 
would be a low density interstellar medium, or a lower mass- 
loss rate during the main sequence. Either would create a lower 
thermal pressure in the main sequence bubble, which in turn 
would place the wind termination shock further away from the 
star. That this can indeed occur would seem to follow from the 
calculations of Chevalier et al. ( 2004 1, which give a wind termi- 
nation shock 28 pc from the central star for GRB 990123. The 
intermediate velocity and 50 kms ' lines can be used to find 
out how much time the star spent as a Wolf-Rayet star before 
producing a supernova. 

The intermediate velocity line is only visible during WRl 
and the first part of WR2, although the absence of this line may 
be the result of the position of the shell fragments (see Sect. 
|5}- The line at 50 kms ' appears later, marking the transition 
from WR2 to WR3, and thus providing another reference 
point for the age of the progenitor star, although we cannot be 



certain of being able to identify this line at all, since it is so 
close to the zero velocity line. 



7.2. The circumstellar medium 

While some gamma-ray burst afterglows indicate a 1/r^ 
density distribution in the circumstellar medium, others appear 
to indicate a constant density (Chevalier et al. l2004t . The only 
areas with constant density in the circumstellar medium are 
the hot bubbles of shocked wind material, which drive the 
shells outward. If the gamma-ray burst afterglow indicates a 
constant density medium, this means that the relativistic jet 
has reached the hot bubble and therefore already passed the 
freely expanding wind region. Consequently, in a constant 
density afterglow profile, the high velocity line caused by the 
Wolf-Rayet wind should not be visible, since the afterglow 
source is already between that area and the observer. Similarly, 
the presence of a high velocity absorption line would seem to 
indicate that the afterglow was caused in a 1 /r^ density region. 
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7.3. GRB 021004 

The best observations so far of the spectrum of a Gamma-ray 
burst afterglow have been made of GRB 021004. Identification 
of the C IV and Si IV lines indicate distinctive absorption lines 
at velocities: kms"\560 kms"', and 3000 kms"' (Schaefer 
et al. I2003> . although the last hne may actually consist of 
two separate components, both at high velocities. This pat- 
tern agrees qualitatively with our column density simulations 
during the early Wolf-Rayet phase, which would indicate that 
the progenitor star only spent a short time as a Wolf-Rayet 
star A faster Wolf-Rayet wind (~ 3000 kms ' rather than 
2200 kms ') remains within acceptable parameters for mas- 
sive stars. Shell velocities of ~ 500 kms"' can be reached by 
adjusting the wind velocities and mass-loss rates of the red su- 
pergiant and Wolf-Rayet winds (cf. Eq. 0}. Calculations done 
by Li & Chevalier ( I2003t indicate that the afterglow of GRB 
021004 results from interaction with stellar wind rather than 
a constant density medium. This agrees with the presence of 
high velocity lines in the spectrum. If there are indeed two high 
velocity lines instead of one, this might be explained by a two- 
stage Wolf-Rayet wind (a fast wind, followed by a slower wind 
as the star expands due to He-shell burning), or by the presence 
of clumps in the Wolf-Rayet wind, or by acceleration of part of 
the Wolf-Rayet wind by radiation from the Gamma-ray burst 
itself. 

VLT spectra of GRB 021004 were presented by Fiore et al. 
<2004> . These observations show a larger number of absorption 
lines at low velocities, leading to a total of six different lines in 
both CIV and Si IV. The large number of intermediate veloc- 
ity components (3) can be explained as the result of multiple 
clumps of the shell fragments appearing in the line of sight (see 
Sect.|5}. The column densities for all these lines are on the or- 
der of 10'"* to 10'^ particles per cm^. In order to compare these 
to our own results, we calculated the total column densities in 
the intermediate velocity and 2200 kms"' lines. The results 
in Fig.ll3lshow that total column densities for the Wolf-Rayet 
wind line are nearly constant throughout the Wolf-Rayet phase. 
Not surprisingly so, since the mass loss rate and wind velocity 
remain unchanged. The location of the wind termination shock 
varies, so the total Wolf-Rayet wind area changes over time, but 
this takes place at a range where the wind density is extremely 
low (see Figs.|3]to|5}, therefore the change in column density 
is minimal. The column density of the intermediate velocity 
feature is strongly time dependent. After all, this feature is the 
result of a moving shell which expands into a spherical volume; 
therefore, the density in the shell will decrease with p ~ r"^ un- 
til the shell remnants hit the outer edge of the bubble. Column 
densities for the lines in Fig^jare on th eorder of 10"^ to 10"^ 
in g cm"^, which corresponds to 10'^ to 10'^ amu per cm"^. As 
carbon has an average weight of about 12 amu, this amounts to 
a particle column density of lO'^ to 10'^ carbon ions per cm"^, 
if the gas were pure carbon. Obviously, this is not the case, so 
the column density of carbon ions would probably be on the or- 
der of 0. 1 to 10 percent, which corresponds quite well with the 
figures found by Fiore et al. ( I2004> . We should keep in mind 
that in Fig. ^] the column densities have been averaged over 
the 200 radial lines. This does not affect the Wolf-Rayet wind 




V, [km s 



Fig. 9. Angle-averaged column density per velocity interval 
dciv,Au) with Av - I kms"' as function of the radial veloc- 
ity at the same moments in time as the circumstellar medium 
plots in Figs.|3]to|5l The formation of the 50 kms"' line can 
be observed, as well as the disappearance of the intermediate 
velocity feature (cf. also Fig. 01. 



line, since the density of the wind is more or less spherically 
symmetric. However, the shell fragments, which make up the 
intermediate velocity line do not have such a symmetry, so col- 
umn densities in the line of sight may vary considerably (see 
Figs.[lO|through[l3. 

A pattern of blue-shifted absorption lines similar to those 
observed in a gamma-ray burst afterglow may be visible in the 
spectra of Type Ib/c supernovae, since these supernovae also 
have Wolf-Rayet stars as progenitors. Dopita et al. ( 1984) at- 
tributed narrow P Cygni profiles seen in Ha and H6 in the 
spectrum of SN 1984E to a Wolf-Rayet wind with a velocity of 
about 3000 kms"'. Similarly, Bowen at al. (2000 1 and Fassia 
et al. (.200 1 ) note the presence of a broad absorption feature 
covering ~ 350 kms"' in the spectrum of SN 1998S, which 
may be caused by a moving shell. 

8. Final remarks 

In this paper, we describe a powerful method for constraining 
the progenitors of GRBs and SNe through narrow circumstel- 
lar spectral lines (see also van Marie et al. 2005 ). While we 
expect a Type Ib/c supernova from practically all massive stars 
whose hydrogen-rich envelope is stripped, only a small fraction 
of these stars — i.e. (0.002 . . . 0.004) (van Putten 2004 1 — pro- 
duces a gamma-ray burst. We have shown that the characteris- 
tics of the circumstellar line spectrum change significantly dur- 
ing the Wolf-Rayet evolution — in line with the fact that only 
a fraction of the Galactic Wolf-Rayet stars are still surrounded 
by a nebula — which means that the circumstellar lines which 
are observable in a GRB afterglow or in a supernova explosion 
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Fig. 10. The presence of gas with a radial column density above 
a threshold value (10"'' g cm"^) moving in three velocity inter- 
vals, as a function of angle and time. The angle is measured 
from the pole down, so an angle of 90 degrees is equivalent to 
the equatorial plane. Red color indicates the presence of gas in 
the 150-300 kms ' interval, while green and blue indicate the 
same for the 400-450 kms"' and 450-600 kms"' intervals, 
respectively. Mixed colors indicate that more than one inter- 
val has a high enough column density; e.g., white indicates the 
presence of gas with high enough column density in all three 
intervals. Pink areas are of particular interest, as these indicate 
a high column density in both the 150-300 km s ' and the 450- 
600 kms"' intervals, whereas the 400-450 kms ' interval has 
a low column density. This means that the intermediate veloc- 
ity feature will show two well-separated components (see also 
Fig.O- 

produced by a Wolf-Rayet star are sensitive indicators of the 
duration of the preceding Wolf-Rayet phase. 

According to stellar evolution models, the duration of the 
Wolf-Rayet phase, which is the final evolutionary phase before 
its death, can vary between almost zero and 10* yr, where it is 
larger for higher mass and/or metallicity (Meynet & Maeder 
12005 1. The fact that only part of the Galactic Wolf-Rayet stars 
actually show a nebula means that those tend to disperse on a 
time scale which is shorter than the mean life time of Wolf- 
Rayet stars, about 5 x lO^yr. Thus, Wolf-Rayet nebulae evolve 
on the evolutionary time scale of their central star. Therefore, 
the corresponding absorption lines at the time the star explodes 
may be used, e.g., to distinguish whether a Wolf-Rayet star died 
young or old. 

In principle, the evolution of main sequence and Wolf- 
Rayet bubbles are influenced — although weakly — by vari- 
ous stellar and interstellar medium parameters (Weaver et al. 
I1977> . However, the life time of the Wolf-Rayet shell, or of 
the intermediate velocity absorption line feature, is on the or- 
der of shell radius divided by shell velocity. For shell velocities 
obtained within our model, this time scale is some lOOOOyr. 
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Fig. 11. The probability of observing the intermediate veloc- 
ity feature as function of time. To estimate this probability, we 
have taken the 2D results during the early Wolf-Rayet phase 
(WRl and the first half of WR2) and checked along each ra- 
dial gridline whether the intermediate velocity line would be 
visible. The criterion for this was a five-fold increase in the 
column density between the area around 100 kms"' and the 
area around 200 km s" ' . 

As the shell velocities coincide well with the values derived 
for the intermediate velocity absorption components in the af- 
terglow spectra of GRB 021004, the conclusion that the Wolf- 
Rayet phase of its progenitor lasted a much shorter time than 
the average Wolf-Rayet life time is in fact independent of the 
particular choice of stellar and interstellar medium parameters 
used in our model. 

A short Wolf-Rayet phase of the progenitor of GRB 021004 
rules out that the underlying Wolf-Rayet star was formed 
through the LBV-scenario (cf. Garcia-Segura et al. I1996at . 
since in this case the duration of the Wolf-Rayet phase cor- 
responds to the full duration of core helium burning (about 
5 X lO^yr), and the nebula would be expected to have dispersed 
by the time the central star dies. Early mass transfer (so called 
Case A or Case B; cf. Wellstein & Langer ll999i . which — just 
like LBV type mass loss — already produces an early Wolf- 
Rayet phase during core helium burning, can also be ruled out. 

The presence of an intermediate velocity component in the 
afterglow spectrum of GRB 021004 indicates that the Wolf- 
Rayet shell was still intact when the star exploded, i.e. that 
the Wolf-Rayet phase was relatively short. This could imply a 
progenitor zero age main sequence (ZAMS) mass close to the 
lower ZAMS mass limit for Wolf-Rayet star formation, which 
is at about 25 Mq for solar metallicity (Meynet & Maeder l2005l 
and larger for smaller Z. Alternatively, the progenitor of GRB 
021004 could have evolved into the Wolf-Rayet phase through 
so-called Case C mass transfer, which occurs at the end of core 
helium burning, and which is assumed to give rise to a common 
envelope phase (cf. Podsiadlowski et al. .l992t . In that case, the 
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Fig. 12. Column density along a single line of sight, i.e. 
dc(y, Ay) with At; = 1 kms"' for an angle of 61° at time t 
= 4.5393 Myr. Angle and time were chosen from Fig. 1101 so 
that two well-separated components of the intermediate veloc- 
ity feature are present. 
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Fig. 13. The total column density of the intermediate velocity 
and 2200 kms ' lines, respectively. The column density of all 
velocity intervals between 150 and 300 kms"' and 2150 and 
2200 kms ' has been added up. Clearly, once the Wolf-Rayet 
wind starts, the total column density within the Wolf-Rayet 
wind remains more or less constant. The intermediate velocity 
feature is extremely time dependent, since this line is caused by 
a moving shell that expands into a spherical volume. Note that 
these column densities where averaged over 200 radial lines. 
This makes little difference for the Wolf-Rayet wind line, since 
the density of the wind is barely angle dependent. The column 
density in the intermediate velocity line, however, is angle de- 
pendent, as is shown in Fig. 1101 



ZAMS mass of the progenitor would have to be smaller than 
the lower ZAMS mass limit for Wolf-Rayet star formation. 

A short duration of the Wolf-Rayet phase might also avoid 
strong angular momentum loss of the stellar core associated 
with the strong Wolf-Rayet mass loss (Langer [l998l Petrovic et 
al. 1200 5V which otherwise might reduce the core angular mo- 
mentum below the critical value required within the collapsar 
model (MacFadyen & Woosley 1999 1. A late common enve- 
lope phase, on the other hand, might even lead to a spin-up of 
the Wolf-Rayet star, either through tidal interaction in a tight 
post-common envelope binary system (e.g. van Putten |2004> 
or through angular momentum gain in the course of a merger 
event (e.g.. Fryer & Heger 2005 1. 

Observed circumstellar spectral lines for a larger number 
of SNe Ib/c or GRB afterglows might enable us to empirically 
constrain the subsample of Wolf-Rayet stars that are capable of 
producing a gamma-ray burst. Since prototype observations ex- 
ist for both cases, we may be able to reach this aim within a cer- 
tain period of time. However, more models of the circumstellar 
medium of massive supernova progenitors at various metallic- 
ities, and for single stars and binary systems, are also needed 
before more general conclusion can be drawn. 
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